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INTRODUCTION 
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An  understanding  of  the  effects  of  environmental  stress  on  human 
performance  is  essential  for  specifying  man-machine  system  design  and 
for  predicting  mission  effectiveness.  Noise,  high  temperature,  accel- 
eration, and  vibration  are  among  the  environmental  stressors  which  affect 
aircrew  members  during  flight.  These  effects  are  seen  as  alteration  of 
physiological  function  and  general  decrements  in  performance  capabilities. 
Laboratory  studies  done  with  man  are  often  compromised  by  the  fact  that 
only  certain  ranges  of  stress  can  be  employed.  One  alternative  is  to 
develop  a series  of  tasks  similar  to  those  of  aircrew  members  or  remote 
pilots  and  train  subhuman  primates  to  perform  such  tasks.  This  alter- 
native offers  the  unique  feature  of  enabling  the  application  of  stresses 
at  all  levels.  It  further  enables  certain  surgical  manipulations 
which  afford  access  to  the  physiological  and  biochemical  factors  of  which 
performance  decrements  under  stress  are  assumed  to  be  a function.  It 
is  only  in  this  way  that  the  fundamental  mechanisms  involved  can  be 
elucidated. 

It  is  also  noteworthy  that  with  the  emergence  of  the  Remote  Pilot 
Vehicle  (RPV)  concept,  attention  becomes  redirected  to  the  tracking  task 
per  se , and  to  the  less  obvious  but  nevertheless  stressful  factors  which 
are  specific  to  the  tracking  task  (e.g.,  fatigue  and  boredom).  Accord- 
ingly, the  laboratory  tracking  situation  has  direct  applicability  to 
the  RPV  paradigm  since  it  closely  resembles  that  situation.  In  fact, 
results  of  tracking  experiments  conducted  in  the  laboratory  can  be 
directly  extrapolated  to  the  RPV  without  the  contaminating  effects  of 
physiological  stresses  due  to  vehicle  dynamics. 
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The  above  mentioned  goals  can  be  achieved  most  efficiently  through 
a viable  program  involving  an  integrated  approach.  One  element  in  such 
an  approach  is  modeling,  a tool  that  can  aid  in  the  identification  of 
basic  mechanisms.  The  results  of  modeling  efforts  are  best  utilized 
when  they  are  closely  tied  to  an  experimental  program.  This  arrangement 
can  produce  an  iterative  process  whereby  the  analytical  results  from 
models  can  serve  as  aids  in  the  interpretation  of  experimental  data, 
guides  in  planning  future  experiments,  and  as  projections  of  results  for 
studies  too  numerous,  costly,  or  difficult  to  conduct.  Similarly,  the 
results  of  the  experimental  program  can  be  used  as  "updates"  for  the 
modeling  studies.  Under  this  ideal  arrangement,  the  iterative  process 
can  proceed  at  a relatively  rapid  rate.  It  is  this  rationale  upon  which 
the  present  research  program  was  developed. 

Specifically,  several  potential  methods  to  produce  tracking  behav- 
ior in  the  Rhesus  monkey  were  explored  in  order  that  the  most  efficient 
procedure  could  be  produced.  Once  developed,  various  parameters  of  the 
method  were  studied  experimentally.  Concomitant  with  this  advancement, 
a mathematical  model  was  developed  and  applied  to  the  data.  A full 
description  of  each  endeavor  appears  in  Sections  1 and  2 respectively 
of  this  report. 


SECTION  1.:  TRACKING  BEHAVIOR 


A.  TRANSITION  FROM  PURSUIT  TRACKING  TO  COMPENSATORY  TRACKING 

In  order  to  appreciate  the  significance  of  this  manipulation,  it 
should  first  be  noted  that  the  early  phases  of  this  project  dealt 
exclusively  with  pursuit  tracking.  This  was  due  to  the  fact  that  the 
original  training  procedures  were  developed  under  another  contract 
(AFOSR  F3361 5-72-C-ll 12)  for  this  task.  Furthermore,  the  pursuit  track- 
ing equipment  had  already  been  constructed,  and  while  the  more  preferred 
compensatory  system  was  under  development,  a series  of  pursuit  tracking 
sessions  was  given  to  three  animals.  The  purpose  of  this  initial  session 
was  to  re-establish  asympototic  performance  levels  which  could  be  used 
as  baseline  data  to  compare  with  those  obtained  in  subsequent  compensa- 
tory tracking  situations. 

All  animals  received  our  standard  series  of  pursuit  tracking  expo- 
sures. Each  daily  session  was  comprised  of  three  15-minute  blocks. 

Each  block  was  sectioned  into  30  second  trials  separated  by  a 30  second 
intertrial  interval.  Although  the  precise  effect  of  this  set  of  time 
parameters  is  unknown,  this  particular  combination  of  time  variables 
has  met  with  considerable  success. 

Prior  to  the  shift,  percent  of  time  on  target  for  the  pursuit  task 
averaged  96.4%,  with  a range  from  92.6%  to  99,8%.  These  values  were 
based  upon  the  last  five  sessions  before  the  shift  in  which  target 
width  was  1.0  inch  with  an  input  signal  of  .122  cps  and  .163  cps  (sum 
of  two  sine  waves). 
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.text,  these  animals  we^e  placed  on  the  compensatory  tracking  task, 
the  basic  parameters  of  which  were  a 1.0  inch  target  and  a random  input 
signal  of  .05  Hz  bandwidth.  Time  on  target  scores  were  collected  for 
each  animal  over  a five  session  block,  pursuit  and  compensatory  scores 
appear  in  the  following  table. 

PERCENT  OF  TIME  ON  TARGET,  TOT  (5  sessions) 

Pursuit  Compensatory 


#065H 

99.8 

94.7 

#065B 

96.9 

99.0 

#476B 

92.6 

69.9 

X 96.4 

X 87.5 

Clearly,  two  of  the  three  animals  manifested  negligible  effects  of  the 
transfer.  In  fact,  one  animal  appeared  to  improve  slightly.  For  the 
third  animal  ( #476B ) the  reduced  time  on  target  score  for  compensatory 
tracking  is  due  solely  to  the  first  session  during  which  this  animal 
"tracked”  at  a mere  54%.  This  monkey  improved  considerably  in  subse- 
quent sessions,  and,  when  the  first  session  was  dropped  from  the 
analysis,  TOT  scores  for  the  animal  rose  to  39.6%.  With  this  corrected 
value  inserted  into  the  analysis,  the  group  mean  was  elevated  to  94.4%, 
which  compares  favorably  with  the  96.4%  value  obtained  for  the  pursuit 
task. 

In  addition,  and  for  comparative  purposes,  two  naive  animals  were 
given  the  standard  response  shaping,  followed  by  a block  of  five  daily 
sessions  on  the  compensatory  task.  Their  time  on  target  scores  ranged 
from  76.4%  to  83.0%  (X  of  87.2%)  for  the  five-day  period. 


The  overall  conclusion  to  be  drawn  from  this  research  is  that  no 
detrimental  effects  are  in  evidence  following  a shift  from  pursuit  to 
compensatory  tracking  when  the  current  set  of  parameters  is  employed. 
This  finding  agrees  with  Poulton's  (1974)  conclusion  that  in  human 
tracking  situations,  the  pursuit-to-compensatorv  shift  can  be  made 
without  decrement.  It  does,  therefore,  illustrate  still  another  simil- 
arity between  the  human  and  subhuman  tracking  task. 

Unfortunately,  the  other  half  of  Poulton's  conclusion  regarding 
asymmetrical  transfer  between  pursuit  and  compensatory  tracking  was 
not  tested  in  this  experiment.  That  is,  because  the  primary  mission  of 
-e  present  research  was  directed  toward  establishing  compensatory 
, acking  in  pursuit  trained  animals,  the  reverse  shift  (compensatorv- 
to-pursuit)  was  not  attempted. 

In  addition,  and  based  upon  the  comparison  with  the  two  control 
subjects,  it  appears  that  not  only  is  the  transfer  from  pursuit-to- 
compensatory  tracking  possible,  the  transfer  is  in  the  positive  direc- 
tion: That  is,  the  three  experimental  animals  appear  to  have  benefited 

from  the  previous  pursuit  tracking  experience  since  their  scores  ranged 
considerably  higher  than  did  those  of  the  control  subjects  which 
received  no  such  training.  Parenthetically,  it  should  be  mentioned  that 
at  the  conclusion  of  this  project,  the  control  animals  were  showing  a 
steady  increase  in  TOT  scores.  Thus,  their  previously  mentioned  scores 
were  not  based  upon  asymptotic  performances.  In  fact,  the  acquisition 
functions  for  these  subjects  are  similar  to  those  of  the  experimental 
animals  on  the  original  pursuit  tracking  task. 

Finally,  the  fact  that  the  control  subjects  learned  the  compensa- 


tory task  at  all,  attests  to  the  fact  that  our  training  procedure  is 
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indeed  applicable  to  either  compensatory  or  pursuit  tracking.  The 
critical  features  of  this  procedure  have  been  spelled  out  elsewhere 
(Lafferty,  Edwards,  McCoy,  McCutcheon,  1973)  and  will  not  be  belabored 
here.  Suffice  it  so  say  that  the  procedure  we  have  developed  appears 
to  fit  either  task  and  that  no  special  retraining  procedures  are  required 
to  train  subhuman  primates  on  either  pursuit  or  compensatory  tracking. 

It  will  be  recalled  that  this  feature,  that  of  a universal  training 
method,  was  one  of  the  primary  goals  of  this  research. 

B.  TITRATED  SHOCK 

Our  standard  training  procedure  employed  continuous  high  levels  of 
shock  at  programmed  rates.  The  presumed  advantage  of  this  technique 
was  that  immediate  feedback  was  provided  when  errors  occurred.  Also, 
the  parameters  involved  (e.g.,  frequency,  intensity)  could  be  easily 
specified.  However,  it  is  also  true  that  shock  administered  in  this 
manner  has  other  properties  which  probably  adversely  affect  tracking 
behavior.  One  obvious  example  is  the  skeletal  reactions  (e.g.,  thrash- 
ing, stick  slamming,  etc.)  which  accompany  shock  and  which  interfere 
with  tracking. 

The  titration  procedure  represented  an  attempt  to  remove  the  above- 
mentioned  emotional  concomitants  which  shock  produces,  and  yet  retain 
shock  as  the  controlling  agent  in  these  studies.  For  present  purposes, 
the  titration  concept  was  applied  as  follows:  When  the  subject  made 

an  error,  shock  was  applied,  but  at  a reduced  (or  perhaps  nonaversive) 
level.  The  intensity  of  the  shock  increased,  however,  so  long  as  the 
error  existed.  When  the  shock  level  became  sufficiently  high,  the  sub- 
ject corrected  the  error  by  returning  the  controlled  element  to  target. 
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In  this  situation,  shock  intensity  was  a direct  function  of  error  time. 
Under  this  condition,  shock  intensity  began  at  1 ma  and  increased  at  a 
rate  of  1 ma  to  a maximum  of  5 ma. 

During  the  latter  phases  of  this  project,  a titration  shocker  was 
purchased,  and  after  several  "false  starts,"  programmed  into  the  com- 
pensatory tracking  system.  One  naive  animal  was  shaped  to  hold,  and 
later  to  move  the  control  stick  in  the  usual  manner.  It  was  then  placed 
on  the  titration  shock  schedule. 

Although  premature  conclusions  are  always  dangerous,  it  does  not 
appear  that  this  subject  was  tracking  as  well  as  some  other  subjects 
which  received  the  standard-shock  treatment.  This  particular  animal 
seemed  to  have  learned  more  slowly,  and  maintained  a consistently  lower 
asympototic  performance  level  than  the  "conventionally  trained"  animals. 
On  the  other  hand,  we  have  learned  that  great  individual  differences 
exist  between  Rhesus  monkey  subjects,  and,  therefore,  another  naive 
animal  was  underqoing  training  according  to  the  same  set  of  parameters 
at  the  termination  of  this  project. 

In  principal,  we  believe  that  the  titration  technique  is  a viable 
one.  At  the  same  time,  it  is  risky  to  draw  conclusions  at  this  point 
in  time.  Additional  subjects  are  required  along  with  certain  schedule 
adjustments  before  firm  conclusions  can  be  drawn.  Unfortunately , funds 
for  tli i s reseach  were  not  available. 

C.  FOOD  MAINTAINED  TRACKING 

Colony  confinement,  handling,  chairing,  confinement  in  the  exper- 
iment chamber  and  shock  are  all  stressful  to  the  Rhesus  monkey.  This 
point  was  alluded  to  in  the  previous  section.  The  overall  effect  of 
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such  factors  may  be  that  a high  stress  level  is  produced  which  interferes 
with,  and  in  some  cases  prohibits,  the  learning  and  maintenance  of  a com- 
plex and  delicate  task  such  as  tracking.  While  the  above-described 
titration  procedure  attempted  to  retain  shock  as  a reinforcing  agent 
through  a reduction  in  shock  intensity,  the  present  section  focuses  upon 
another  source  of  control,  food. 

The  basic  features  of  this  procedure  involved  training  the  subject 
to  maintain  a controlled  element  on  target  for  progressively  longer 
time  periods  in  order  to  obtain  food  pellets.  I:  order  to  avoid  satia- 
tion, the  subject  was  given  food  pellets  only  at  programmed  intervals 
during  the  tracking  session.  At  other  times,  a tone  was  substituted 
for  food.  By  virtue  of  its  previous  association  with  food,  the  tone 
has  acquired  the  capactiy  to  reinforce  behavior  secondarily.  Actual 
tone-food  pairing  occurred  on  a fixed  ratio  basis  of  5:1.  For  each 
criterion  time  on  target,  the  tone  was  presented.  Food  accompanied 
the  tone  (and  maintained  its  strength  as  a conditioned  reinforcer) 
on  every  fifth  presentation. 

The  first  animal  studied  was  trained  jointly  on  food  and  shock. 

When  these  reinforcing  agents  were  later  made  independent,  it  became 
clear  that  shock  was  the  overrriding  agent:  That  is;  the  subject  would 
not  work  for  food  independently  of  shock,  but  would  track  to  avoid 
shock  independently  of  food. 

Later  in  the  grant  period,  a second  animal  was  pretrained  on  the 
FR5  tone-food  pairing  and  subsequently  placed  in  the  compensatory 
trackinn  situation.  After  a few  sessions  this  animal's  behavior  became 
highly  erratic,  and  TOT  scores  diminished  rapidly.  Mild  food  depriva- 
tion did  not  remedy  the  situation.  It  appeared  that  the  conditioned 
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(tone)  was  losing  its  reinforcing  effectiveness.  At  this  tine,  tracking 
^ sessions  were  halted  and  tone-food  pairings  were  resumed.  In  spite  of 

this  additional  training,  the  tone  did  not  appear  to  be  an  effective 
reinforcer.  It  was  also  discovered  that  the  animal  was  consuming  pro- 
gressively fewer  sucrose  pellets.  Finally,  it  would  not  longer  eat  them 
when  made  freely  available  in  the  home  cage. 

Other  animals  have  evidenced  similar  behavior  with  other  sorts  of 
reinforcing  agents.  The  decline  in  the  reinforcing  effectiveness  of 
these  conventional  reinforcers  is  undoubtedly  responsible  for  the  highly 
erratic  data  yielded  in  the  food-reinforcement  situation. 

The  fact  that  we  have  not,  at  this  point,  established  clearly  that 
compensatory  tracking  can  be  maintained  on  a food  reinforcement  schedule, 
does  not,  in  our  opinion,  vitiate  against  the  basic  premise  of  the  pro- 
cedure. Any  behavior  is  only  as  strong  and  reliable  as  are  the  conditions 
used  to  support  that  behavior.  If  the  reinforcers  upon  which  the  behav- 
ior depends  lose  their  effectiveness,  the  behavioral  output  will  dimin- 
ish. This  does  not,  however,  reduce  the  utility  of  the  food-tracking 
premise.  Rather,  it  simply  indicates  that  other  reinforcing  agents  must 
be  employed.  At  the  conclusion  of  the  grant  period,  we  were  of  the 
opinion  that  this  problem  had  been  solved  but,  due  to  a lack  of  support- 
ing funds,  we  were  not  able  to  pursue  the  problem  further. 

D.  HUMAN  VS.  SUBHUMAN  OPERATORS  IN  TWO  TRACKING  SITUATIONS 

As  part  of  our  program  to  study  the  comparative  similarities  and/or 
differences  between  human  and  monkey  operators  under  identical  situations, 
data  were  collected  from  these  two  species  for  both  pursuit  and  compen- 
satory tracking  tasks.  In  this  connection,  it  is  of  interest  that  a 
paper  appeared  recently  (Bachman,  Jaeger,  and  Newsom  1976)  which  focused 
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upon  the  issue  of  man-monkey  comparisons  in  a tracking  situation.  Using 
an  argument  similar  to  that  which  we  have  employed  for  some  time  (i.e. 
that  a monkey  model  is  necessary  for  the  application  of  stress  to  man), 
these  authors  argued  convincingly  that  we  must  first  demonstrate  that 
these  species  do  show  similar  performances.  Based  upon  a primate-train- 
ing task  quite  similar  to  that  discovered  in  our  laboratory  some  four 
years  ago,  Bachman,  et  al . produced  results  which  do  lend  credibility 
to  the  monkey-man  extrapolation  since  the  two  species  did  perform  simi- 
larly on  the  task  studies.  However,  their  results  were  based  solely 
upon  the  compensatory  tracking  situation.  Our  present  results  both 
support  and  extend  those  of  Bachman  et  al . , since  our  data  are  based 
upon  both  compensatory  and  pursuit  tracking  situations. 

In  the  initial  phase  of  the  investigation,  four  humans  and  three 
Rhesus  monkeys  were  studied  in  the  pursuit  tracking  situation.  Target 
size  was  1.0  inches  and  a sum  of  two  sine  waves  (.122  and  .163  cps) 
was  used  as  the  input.  Trials  were  30  seconds  in  duration  and  separated 
by  a 30  second  intertrial  interval.  Three  fifteen-trial  blocks  were 
given  to  each  subject  in  a single  session.  For  the  monkeys  "off  target' 
performances  were  punished  with  shock.  For  the  humans,  these  errors  were 
punished  on  an  error-cost  basis  which  was  subtracted  from  the  total  amount 
of  money  ($3  maximum)  which  they  could  earn.  Next,  all  subjects  were 
studied  in  the  compensatory  tracking  situation,  characterized  by  a target 
size  of  1.0  inches  and  an  input  bandwidth  varying  from  .05  to  .15  cps. 
Error  contingencies,  trial  durations  and  intertrial  intervals  were  as 
in  the  former  task. 

Although  the  fine-grained  behavioral  outputs  remain  to  be  analyzed 
and  mathematized,  a preliminary  evaluation  can  be  made  from  the  following 
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table  based  on  percent  of  tirne-on-target  scores  (TOT). 


1 

Pursuit 

Compensatory 

.05 

.15 

Humans 

99.8 

99.9 

99.3 

I 

99.6 

99.5 

99.3 

1 

98.5 

98.9 

98.5 

96.7 

97.6 

97.4 

Monkeys 

99.3 

94.7 

94.2 

96.9 

99.0 

98.4 

92.6 

89.6 

90.8 

From  the  above  table,  it  is  clear  that  both  human  and  subhuman  op- 
erators perform  quite  well  on  both  tasks.  Very  few  errors  occur  in  any 
situation.  More  importantly,  both  humans  and  monkeys  show  negligible 
effects  of  the  parameter  changes.  The  functional  relationship  between 
TOT  and  parameter  changes  is  quite  similar  for  each  individual  organism 
studied.  It  is  the  point  which  is  of  critical  significance;  viz  that 
the  functional  relationships  are  the  same,  regardless  of  the  species . 

The  results,  along  with  those  of  Bachman  et  al . , offer  strong 
support  for  the  contention  that  human  and  monkey  operators  are  similar 
on  the  two  tasks.  Any  differences  which  exist  are  quantitative,  not 
qualitative.  Therefore,  the  emotionally-based  argument  offered  by  some 
that  man  and  monkey  are  functionally  different  is  not  supported  by 
empirical  fact. 
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E.  EFFECTS  OF  TRANQU ILIZING  DRUGS  OR  PURSUIT  TRACKING 

Before  discarding  the  mechanical  pursuit  tracking  facility  in  favor 
of  the  compensatory  system,  two  animals  were  given  a series  of  tracking 
sessions  under  various  dosages  of  two  tranqui 1 izing  agents,  chlorpro- 
mazene  (CPZ)  and  pentobarbital.  CP7.  is  classified  as  an  antipsychotic 
drug  and  is  viewed  as  a "major  tranquilizer".  Pentobarbital  is  a barbi- 
tuate,  and  is  classified  as  a "minor  tranquilizer".  Relatively  little 
is  known  regarding  the  effectiveness  of  these  drugs  on  complex  pyscho- 
motor  activities  such  as  tracking.  Therefore,  one  purpose  of  the  study 
was  to  determine  whether  the  two  chemical  agents  would  have  different 
effects  on  tracking  ability.  Another  purpose  was  related  to  the  earlier- 
mentioned  "emotional  effects"  of  high  stress  levels  associated  with  the 
shock-maintained  tracking  situation.  Specifically,  it  seemed  possible 
that  low  dosages  of  tranquilizers  could  actually  improve  tracking 
efficiency. 

As  with  many  psychopharmacol ogical  experiments,  the  general  plan 
was  to  establish  a baseline,  administer  a drug  at  a low  dose  level, 
and  gradually  increase  the  dose  on  subsequent  drug  sessions  until  the 
performance  level  changed  significantly.  Next,  for  purposes  of  repli- 
cation, the  lower  case  levels  we^'o  repeated.  In  the  present  situation, 
drugs  were  given  intramuscularly  at  irregular  intervals,  but  never  more 
than  one  each  week.  Several  non-drug,  or  placebo  sessions  were  always 
spaced  between  any  two  drug  conditions.  One  animal  (Hoppy)  was  given 
CPZ  at  several  dosages  and  the  pentobarbital  at  several  dosages.  The 
drug  order  was  reversed  for  the  second  animal,  Big  Boy.  Both  subjects 
received  the  "standard"  tracking  conditions  described  earlier  (i.e., 


30  second  trials,  30  second  IT  I , ful 1 -intensi ty  shock). 


Computer  analysis  and  mathematical  description  of  the  "fine  grain" 
aspects  of  the  drug-produced  behaviors  is  yet  to  be  accomplished.  Never- 
theless, time  on  target  scores,  TOT  ( in  terms  of  percentages),  appear 
below.  It  is  important  to  realize  that  the  drug  dosages  were  given  in 
the  order  in  which  they  appear. 

CPZ  (Hoppy)  BASELINE  = 99.3% 

Dose  (Mg/ Kg ) .2  .5  1.0  .75  .5  .35  .2 

% TOT  97.1  97.7  41.2  32.8  41.0  54.9  51.8 


CPZ  (Big  Boy)  BASELINE  = 96.3°/ 

Dose  (Mg/Kg)  .2  .35  .2  .15 

% TOT  92.2  40.0  76.2  95.2 


PENTOBARBITAL  (Hoppy)  BASELINE  = 99.811 
Dose  (Mg/Kg)  5.0  8.0  1.20  3.5  5.0 

% TOT  94.0  74.7  71.4  58.1  59.3 


_ PENTOBARBITAL  (Big  Boy)  BASELINE  = 96.3% 

Dose  (Mg/Kg)  5.0  8.0  12.0  15.0  8.0  5.0 

% TOT  90.5  88.0  90.0  41.0  86.3  89.1 

Looking  first  at  the  effect  of  CPZ,  the  lowest  dose  (.2  Mg/Kg)  had 
only  a slight  effect  on  tracking  efficiency.  This  was  true  for  both 
animals.  At  this  point,  the  similarity  of  their  reactions  to  the  drug 
ended.  For  one  animal  (Hoppy)  increased  doses  up  to  1.0  Mg/Kg  produced 
systematic  decreases  in  TOT  scores.  More  interesting,  however,  is  the 
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fact  that  the  return  to  lower  dose  levels  was  not  accompanied  by  increases 
in  tracking  efficiency.  Thus,  this  animal's  response  was  dose  dependent 
until  tracking  had  become  severly  attenuated,  but  it  was  dose  independ- 
ent beyond  that  point. 

A different  picture  emerged  with  the  second  animal  (Big  Boy) 
studied  under  CPZ.  In  the  first  place,  this  monkey  was  clearly  more 
sensitive  to  the  drug.  A dose  of  .35  Mg/Kg  was  sufficient  to  suppress 
his  tracking  efficiency  below  the  chance  level.  More  interestingly, 
this  animal's  behavior  appears  to  be  more  dose-dependent.  The  return 
to  lower  doses  yielded  progressive  increases  in  tracking  efficiency. 

The  only  discrepency  in  this  generalization  occurred  at  .2  Mg/Kg  where 
TOT  was  16%  less  on  the  second  administration  of  this  dose.  This  was 
due  mainly  to  a single  episode  in  which  the  animal  removed  his  hands 
from  the  stick. 

Consider  now,  the  effect  of  pentobarbital.  Clearly,  much  larger 
doses  of  this  agent  were  required  to  affect  behavior  measurably.  For 
one  animal,  Hoppy,  the  effect  was  much  the  same  as  with  CPZ;  that  is, 

TOT  scores  decreased  progressively  with  increased  dose  level.  Beyond 
this  point,  however,  decreases  in  dosage  were  not  accompanied  by 
increases  in  tracking  efficiency.  TOT  scores  at  3.5  and  5.0  Mg/Kg  were 
both  considerably  below  that  at  the  initial  5.0  dose  level.  Again, 
this  animal's  tracking  efficiency  was  at  first  decreased  by  increasing 
dosages,  but  later,  TOT  remained  low  and  was  unrelated  to  dose  magnitude. 
In  sharp  contrast,  was  the  effect  of  pentobarbital  on  the  second  animal, 
Big  Boy.  Again,  the  effect  of  dose  strength  on  tracking  efficiency  was 
inverse;  as  dose  level  increased,  tracking  efficiency  decreased. 

Clearly,  the  effect  of  both  drugs  was  dependent  as  much  on  the 
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individual  animal  as  it  was  on  the  specific  dose  employed.  The  most 
interesting  case  is  Hoppy  who  was  unable  to  recover  his  tracking 
efficiency  at  any  dose  level.  Still,  this  subject  did  track  effective- 
ly on  non-drug  days.  As  a matter  of  fact,  this  animal's  performance 
actually  improved  slightly  over  successive  baseline  sessions.  Perhaps, 
for  this  animal,  any  quantity  of  drug  functioned  as  a discriminative 
stimulus  and  its  performance  became  "state  dependent"  (e.g.,  Overton, 
1967).  This  finding  may  also  reflect  the  possibility  that  the  two 
animals  studied  used  quite  different  strategies  in  learning  and  perform- 
ing in  the  tracking  situation,  and  that  the  drugs  studied  here  select- 
ively affected  one  of  these  mechanisms. 
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SECTION  II 

DATA  ACQUISITION,  ANALYSIS  AND  MATHEMATICAL  MODELING 

This  section  describes  the  development  of  a 5 parameter  model  to 
fit  the  primate  data  in  a compensatory  tracking  task.  A computer  iter- 
ative scheme  of  determining  the  optimal  parameters  of  a given  model 
based  on  experimental  data  is  also  detailed.  This  section  begins  with 
a description  of  the  electronic  system  used  in  tracking  experiments 
followed  by  a overview  of  process  control,  data  acquisition  and  analysis, 
and  model  identificaiton  scheme  and  modeling  results. 

A.  TRACKING  SYSTEM 

The  electronic  system  used  in  tracking  experiments  was  custom- 
designed  and  built  at  Wenner-Gren  Research  Laboratory.  The  machine  is 
fully  programmable  in  that  all  parameters  like  "Test",  "Rest"  and  "Grace" 
periods,  and  number  of  "Test/Rest"  cycles  can  be  externally  prepro- 
grammed. In  addition,  the  machine  can  be  used  either  in  the  Pursuit  or 
Compensatory  mode. 

The  flow  chart  of  Fig.  2.1  gives  a block  diagram  description  of  the 
system.  The  system  primarily  consists  of  hybrid-analog  and  digital  sig- 
nal conditioners  designed  to  control  the  X,  Y and  Z axis  display  of  a 
Cathode  Ray  Tube  (CRT).  The  stationary  pattern  of  the  CRT  screen,  shown 
in  Fig.  2.2a,  is  generated  internally  using  analog  circuits.  The  scheme 
used  to  do  this  is  shown  in  Fig.  2.2b. 

Notice  that  the  pattern  of  Fig.  2.2b  can  be  generated  on  a CRT  by 
simply  using  r-j  Sin  wt  and  r^  Sin  wt  on  horizontal  axis  display  control 
and  r^  Cos  wt  and  r^  Cos  wt  on  vertical -axi  s display  control.  Both  r-j 
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Sin  ut  and  r ^ Sin  ut,  and  r-|  Cos  wt  and  r9  Cos  ut  are  tiine-mul tiplexed 
using  Z-axis  control.  However,  if  r9  on  the  vertical  axis  is  greatly 
amplified  compared  to  r ^ on  the  horizontal  axis,  the  pattern  of  Fig. 

2.2a  would  be  generated.  Depending  on  the  mode,  pursuit  or  compensatory, 
the  analog  input  signal,  stimulus  x(t),  is  used  to  control  the  motion 
of  vertical  lines  or  the  central  cursor  (dot)  along  the  X-axis,  respect- 
ively. The  stick  output  y(t)  controls  the  motion  of  the  cursor  both  in 
pursuit  and  compensatory  modes.  An  error  signal  e(t)  is  generated  when 
the  cursor  crosses  boundary  limits  which  turns  shock  on  and  simultaneously 
causes  the  cursor  to  blink.  All  these  signals  x(t),  y(t)  and  e(t)  are 
internally  conditioned  to  lie  within  0-2  volts  and  are  made  available  to 
the  Raytheon-704  computer  for  digitization. 

The  digital  section  is  basically  a quad-time-state  logic  machine. 
External  manual  controls  "Clear"  and  "Go"  enable  the  machine  to  be  reset, 
and  initialized  before  starting  the  experimental  session.  "Clear"  sig- 
nals the  machine  into  the  "Ready"  state  and  simultaneously  triggers  the 
CRT  pattern  of  Fig.  2.2a.  The  "Go"  signal  takes  the  machine  to  "Test" 
state.  In  this  state  the  machine  waits  for  "Grace"  period  T-j  , and  then 
checks  to  see  if  the  subject  is  on  target.  In  an  off  target  condition, 
the  cursor  is  caused  to  flash  and  shock  is  applied  to  the  subject. 

Shock  continues  until  the  subject  is  back  on  target  at  which  time  both 
light-blinks  and  shock  signals  are  turned  off.  During  the  entire  "Test" 
period  the  machine  keeps  track  of  time  and  at  the  end  of  T^  seconds 
"Rest"  state  is  entered  and  all  activity  ceases  for  T^  seconds.  The 
"Rest"  state  is  followed  by  "Increment"  state  in  which  the  machine  checks 
if  preset  numbers  of  "Test/Rest"  cycles  have  been  completed,  if  not, 
a new  cycle  is  initiated.  At  the  end  of  the  last  cycle  the  machine  halts. 
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The  “Grace"  period  T-j  may  be  preset  as  a multiple  of  1 second. 
“Test"  and  "Rest"  periods  T^  and  T^  respectively  can  be  preprogrammed 
as  multiples  of  10  seconds. 

1.  Random  Inputs:  The  input  used  for  exciting  the  system  (causing 

the  cursor  to  move)  was  Gaussian-randoni  in  nature  and  was  derived  from 
a Random  Signal  generator  (Hewlett  Packard).  The  primary  reason  for 
using  random  input  signals  was  to  remove  the  element  of  predictability 
from  the  operator's  response.  In  this  research,  the  analog  random  sig- 
nals of  frequency  bandwidths  0.05  Hz  and  0.15  Hz  were  used. 

B.  EXPERIMENT 

The  main  body  of  the  experiment  consisted  of  a series  of  runs  cov- 
ering the  pursuit  task  initially  and  then  gradually  changing  to  compen- 
satory tasks.  The  range  of  input  frequencies  was  very  low,  only  extend- 
ing to  0.15  Hz.  This  was  done  to  establish  a data  base  in  the  lowest 
frequency  range,  disregarded  thus  far. 

Both  in  pursuit  and  compensatory  tracking  tasks  the  animals  were 
restrained  in  a chair  in  full  view  of  a control  stick  and  CRT  target 
di splay. 

1.  Pursuit  Tracking:  The  two  vertical  lines  of  Fig.  2.2a  were 

moved  either  in  a regular  or  random  fashion,  along  the  X-axis  under  the 
influence  of  an  external  stimulus.  Test  animals  were  required  to  maneu- 
ver the  control  stick  to  keep  the  cursor  between  the  two  vertical  lines. 

2.  Compensatroy  Tracking:  The  external  stimulus  was  applied  to 

the  cursor  forcing  it  outside  the  vertical  lines.  The  test  animals  were 
required  to  manipulate  the  stick  to  position  the  cursor  within  the  ver- 
tical lines. 
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C.  DATA  ACQUISITION  AND  ANALYSIS 

During  each  experimental  session,  random  input  x(t)  to  the  con- 
trolled element  and  stick  output  y(t)  were  digitized  every  16  msecs  on 
the  Raytheon  704  digital  computer.  The  computer  was  programmed  to  re- 
cognize "Rest"  periods  and  discontinue  data  acquisition.  In  brief,  data 
was  recorded  for  15  minutes  on  magnetic  tapes  for  off-line  analysis  on 
the  IBM  370/165  computer.  The  sampling  rate  of  16  msec  was  chosen  to 
include  all  higher  frequencies  in  the  subject's  response. 

For  data  analysis  purposes,  the  magnetic  tapes  were  taken  to  the 
Computer  Center  for  evaluating  a describing  function  using  programs 
based  on  Time-series-analysis  (Gupta,  1977).  This  program  was  developed 
for  use  on  IBM  370/165,  using  subroutines  from  literature  (Reid,  1969). 

The  first  part  of  the  program  combines  subroutines  FRTRM,  SORTG, 
CRSCOR  and  PSD,  to  evaluate  cross-spectral  estimates  (cross-spectra 
between  input  and  output)  and  ; ^ e (cross-spectra  between  input  and  error) 

$.  and  $.  estimates  were  used  to  calculate  the  describing  function  as 
to  le  3 

given  by 

- (jw) 


2.1 


Y(jcu)  = 


10 


*ie  (J*>) 


Since  in  Equation  2.1,  both  and  are  complex  numbers , they  can 
be  combined  to  express  Y(jio)  as  a series  of  real  and  imaginary  parts. 
Re  [*.  ] + j Im  [>. J 


Y(ju>)  = 


lo- 


io- 


Re  Oi0]  + j Im  01e] 

Re  OiQ]  + j Im  [tiQ]  Re  [*1e]  - j Im  [>ie] 


Re  [*1e]  + j Im  [*1fi] 


Re 


Oie]  - J Im  [fie] 


2.2 
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(Re  [>io])2  + (Im  Oie])2  + j(Re[*ie]  . Im[>io] 

" -Re[»i0].Im[»ie]) 

(Re  [*ie])2  + (Im  [^ie])2 

= Y1  (u)  + j Y"  (ui) 

D.  MODEL  IDENTIFICATION  SCHEME 

The  last  phase  of  the  analysis  was  matching  the  describing  function 
estimates  with  a suitable  model.  Model  identification  starts  with  the 
assumption  of  a known  model  based  on  a visual  judgement.  Since  the  model 
transfer  function  G(jw)  can  be  expressed  as  a complex  number: 

2.3  G(ju>)  = G'(u>)  + j G"U) 

an  alogrithm  was  devised  to  determine  optimal  values  of  various  model 
parameters  to  minimize  the  sum  squared  error  S defined  by 

2.4  S (A)  = E (G!  - Y!)2  + (GY  - YY)2 


where  i is  a frequency  index,  G.!  and  GY  are  data  computed  for  the 
assumed  model  and  Yj  and  YY  are  experimentally  observed  data.  G.I  and 
GY  are  computed  for  the  chosen  model  using  arbitrarily  chosen  values  of 
unknown  parameters  A,  the  dimension  of  which  depends  on  the  number  of 
parameters  in  the  model. 

Differentiation  of  Equation  2.4  with  respect  to  a given  unknown 

parameter  A.  is  given  by 
J 


2.5 


2 E (G! 


+ 2 (GY 


YY) 


aGY 


Since  S is  a function  of  A,  a Taylor  series  expansion  in  terms  of  known 
values  of  S and  may  be  made.  (A  is  the  initially  chosen  parameter  set) 
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2.6  S (A)  = S (A°)  + ~ t|—  Aj  + higher  order  trms 

J 

Thus  if  higher  order  terms  in  Equation  2.6  can  be  neglected,  an  inter- 
active program  using  Equations  2.2  and  2.3  may  be  devised  such  that  S(A) 
tends  to  zero  (perfect  fit).  The  program  begins  with  an  arbitrarily 
chosen  parameter  vector,  A°.  The  upper  and  lower  limits  on  the  unkonwn 
parameters  are  required  both  to  limit  the  search  area  and  to  physically 

i C 

interpret  the  data.  The  function  value,  S(A°),  G , G"  and  are  com- 

1 1 ' t\  • 

0 

puted  using  Equation  2.6  and  the  chosen  model.  If  S(A°)  and  the  magni- 
tude of  the  gradient,  G(A°),  defined  by, 

2-7  = <ff°>  ■ $7  - % • )T 

are  less  than  a specified  tolerance  the  computer  search  for  a local 
minima  is  over  and  a new  one  can  be  made  with  a different  starting 
parameter  vector.  This  is  usually  not  the  case  whence  a move  in  the 
negative  gradient  direction  is  made  to  define  a new  parameter  vector 

2.8  A11  = A ° - a G(A°) 

where  a is  defined  by 


< G(A° ) , G(A°)  > 

Equations  2.8  and  2.9  define  the  steepest  descent  formula  (McGhee,  1967). 
At  the  new  vector  A11  the  curve  fit  error  S(An)  is  computed  and  compared 
with  S(A°).  If  the  error  has  been  reduced  the  new  values  are  assumed 
and  Equations  2.8  and  2.9  are  used  recursively.  In  the  event  the  steep- 
est descent  method  fails,  a second  order  gradient  technique  known  as  the 
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Newton-Raphson  technique  (McGhee,  1967)  has  been  incorporated  which  has 
good  convergence  properties  in  the  vicinity  of  a minima.  The  new  vector 
A11  with  this  formula  is  given  by 

2.10  A11  = A ° - 0T  $ + 0T  o]'1  G(A°) 

where  $ and  a are  rectangular  matrices  of  partial  derivatives  of  error 
terms  (G.!  - Y.!)  and  (GY  - YY)  with  respect  to  the  unknown  parameters. 

The  process  of  iterative  minimization  of  sum-squared  error  continues 
until  no  further  reduction  in  S is  possible  with  both  the  above  gradient 
techniques . 

E.  VERIFICATION  OF  DATA  REDUCTION  TECHNIQUES 

In  order  to  test  the  program's  ability  to  measure  a linear  system 
transfer  function,  three  types  of  simulations  were  made.  The  first  sim- 
ulation to  test  auto-power-spectral  estimates  was  effected  by  digitizing 
two  sinusoidal  waveforms  of  frequencies  2 Hz  and  6 Hz.  The  program  iden- 
tified the  predominant  frequencies  close  to  2 Hz  and  6 Hz.  Fig.  2.3 
shows  power  spectral  estimates  predicted  by  the  program. 

The  second  simulation  of  Fig.  2.4  was  aimed  at  identifying  an  open 
loop  system.  The  open  loop  system  was  simulated  digitally  using  both 
the  standard  z-transform  and  bilinear  z-form  techniques.  The  z-trans- 
form  equation  for  the  system  is  given  by 

o n Y(z)  _ -10  _ -1 0T 

*tTT  = s-a  z = esT  - 1.e-aTz.1 

where  a = 3 

T = 16  msec 


The  bilinear  z-form  equations  for  the  system  of  Fig.  2.3  are 


7.0E+6 
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Fiq.  2.3  Frequency  Identification 


Cy^-ferv,  (closed  Uop)  Xdenilficah** 
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2.12 

and 

2.13 


S = 


1-2 

1+2 


-1 

-7 


2 

T 


V = tan 
2 1 


Actual  and  observed  magnitude  and  phase  plots  are  shown  in  Figs. 

2.6a  and  2.6b.  A comparison  of  the  results  was  made  with  a program 
developed  at  the  University  of  California  (BMD  02T).  It  was  observed 
that  our  programs  fitted  the  actual  results  to  within  an  accuracy  of 

2 

2%.  The  results  are  shown  in  Fig.  2.7  where  the  correlation  index  p 

is  computed  as  a function  of  frequency.  It  should  also  be  noted  that 
2 

the  cumulative  p is  very  nearly  .99,  indicating  that  the  system  is 
linear.  In  order  to  verify  the  program's  ability  to  identify  an  unknown 
system  in  a closed  loop  configuration , the  system  shown  in  Fig.  2.5  was 
digitally  simulated  using  Equations  2.14  and  2.15  to  generate  Y(kT)  and 
E(kT)  respectively. 

2.14  Y(kT)  = {1.323  x(k-2)  - x(k)  + 7.92410S8  Y(k-l) 

-(4.0379456  -.1232coD)  Y ( k-2 ) } / (4.0379456  + .1232^) 

2.15  E (kT)  = { (4.1  + .1232^)  x(k)  - 8.0  x(k-l)  + (4.0  - .1232^) 

x(k-2)  + 7.9241088  e(k-l)  - (4.0379456  - .1232^) 

Y (k-2)  } / (4.0379456  + .1232^) 

Where  is  the  digital  equivalent  of  the  analog  cutoff  frequency  defined 
by  Equation  2.13.  The  results  are  shown  in  Figs.  2.8a  and  2.8b. 

The  program  was  found  convergent  for  various  starting  values  chosen 
for  a known  4-parameter  model  of  Equation  2.16. 

e"TS  (1  + T,s) 

2.16  G(ju)  — 

(1+T2s)(l+T  s) 


3UlllINLJVh 


Fin.  2.6a  Linear  System  (Open  Lood)  Identification 
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The  model  parameter  values  were  chosen  as  x = 0.2  sec,  T-j  = 0.2  sec, 

><  T,  = 0.5  sec,  and  = 5.0  sec.  The  program  identified  these  parameters 

as  r = 0.19997,  T-|  = 0.19997,  T^  = 0.50005  and  T^  = 4.9950  starting  with 
the  arbitrarily  chosen  values  of  x = 0.5,  T-j  = .5,  T^  = 2.5  and  T^  = 5.0. 
The  model  identification  scheme  was  also  used  for  a 5-parameter  model 
of  Equation  2.17.  The. reason  for  choosing  the  model  of  Equation  2.17  is 
the  general  acceptance  of  this  model  in  literature  for  human  subjects. 

Ke-TS  (1  + T,s) 

2.17  G(ju)  = — — - ' 

(1  + T s)(l  + T s) 


With  parameters  chosen  as 

K = 2.0 
x = 0.2  sec 
T^  = 0.2  sec 
T2  = 0.5  sec 
T3  = 5.0  sec 

A typical  set  of  parameters  identified  was 

x = 0.20012  sec 
K = 1 .9990 
T-j  = 0.2005  sec 
T2  = 0.5008  sec 
T3  = 4.995  sec 

_5 

with  s = 0.31970  x 10  . Several  starting  parameter  vectors  were  chosen 

for  both  4 parameter  and  5 parameter  models  and  results  were  comparable 
to  those  discussed  above.  A typical  run  with  6 trials  took  60  seconds 
of  computer  time  i.e.,  10  seconds  per  trial. 

Table  2.1  and  2.2  contain  the  initial  and  final  vectors  for  4-para- 


meter  and  5-parameter  models  respecti vely . 
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Table  2.1  - 4-Parameter  Model 


# of  Initial  Final  Correlation 

Run  $ Iterations  Vector  Vector  Error 
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Table  2.2  - 5-Parameter  Model 


Run  # 

# of 

Iterations 

Initial 

Vector 

Final 

Vector 

Correlation 

Error 

1 

140 

T = . 5 

T “ 

.2001 

.1349  x 10"b 

T!  = -5 

T1  = 

.2005 

T2  = 2.5 

T2  = 

.5003 

T3  - 5.0 

T3  = 

4.995 

K = 1.0 

K = 

1 .9999 

2 

104 

T = .3 

T = 

.1998 

.1445  x 10'6 

T1  = -4 

T1  = 

.1994 

T2  = 2.0 

To  = 

C 

.4991 

T3=6-° 

T3 

5.0041 

K = 2.0 

K = 

2.0001 

3 

42 

T = .6 

T ” 

.2001 

.1491  x 10~6 

T1 = -6 

T1  = 

.2005 

T9  = 4.0 

T2  = 

.5009 

T3  = 6-° 

T3 

4.9957 
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F.  MODELING  RESULTS 

In  the  investigation  detailed  here,  three  Rhesus  monkeys.  Butch, 
Hoppy  and  Big  Boy  were  used  as  subhuman  operators  in  a shock  controlled 
compensatory  tracking  situation.  For  each  of  the  animals,  five  taDes 
were  randomly  chosen,  both  in  .15  Hz  and  .05  Hz  ranges.  One  of  the 
three  monkeys,  namely  Big  Boy,  was  trained  for  only  a few  days  on  0.05 
Hz.  Consequently,  the  number  of  recordings  made  were  not  enough  to  es- 
timate his  describing  function  in  the  .05  Hz  range. 

The  describing  function  data,  namely  magnitude  in  dB  and  phase  angle 
in  degrees,  for  the  three  primates,  is  presented  in  Tables  Ml  through 
M5.  Tables  M2,  M3  and  M4  contain  estimates  for  frequency  bandwidth  of 
.15  Hz  and  Tables  Ml  and  M5  for  .05  Hz.  The  last  two  columns  in  each 
of  the  Tables  M1-M5,  both  for  magnitude  and  phase  angle,  contain  the 
mean  for  the  five  sets  of  data  together  with  the  standard  deviation. 

Both  for  the  .15  Hz  and  0.05  Hz  signal  bandwidths,  magnitude  of  the 
describing  function  averaged  between  -12.0  dB  and  30.0  dB  whereas  the 
phase  angle  averaged  between  10°  to  -115°.  In  the  frequency  range  of 
.15  Hz,  the  standard  deviation  in  magnitude  and  phase  angle  were  re- 
spectively of  the  order  of  0-3  dB  and  3°-30°.  In  the  0.05  Hz  bandwidth, 
however,  the  deviation  was  relatively  higher.  This  was  mainly  due  to 
the  aliasing  introduced  by  the  very  low  frequency  components  in  the 
input  signal.  Mean  describing  function  estimates  together  with  the 
standard  deviation  are  plotted  in  Figs.  2.9  through  2.13. 

Similarly,  Tables  HI  and  H2  contain  the  describing  function  data 
for  human  subjects  Human  2 and  Human  3 for  both  frequency  ranges  0.05 
Hz  and  0.15  Hz.  The  magnitude  and  phase  angle  for  the  human  subjects 
ranged  from  5 dB  to  45dB,  and  -10°  to  -110°  respectively.  It  should 
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PHRSE (DEGREES]  *10'  MAGNITUDE  (DB) 
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Fig.  2.11 
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DESCRIBING  FUNCTION  - BIGBOY 
SIGNAL  BANDWIDTH  --  .15  HZ 

$ - OBSERVED 


I 


Fig.  2.13 
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Table  H 1 - Describing  Function  Data  - Human  Subjects 
Frequency  Range:  0.15  Hz 


Subject 

- H 3 

Subject 

- H 2 

Freq.  (Hz) 
★ 

245.6 

flag.  (dB) 

Phase 

(Degrees) 

Mag.  (dB) 

Phase 

(Degrees) 

1 

27.19 

-62.25 

23.13 

-10.57 

3 

21  .95 

-62.39 

19.10 

-40.38 

5 

21.57 

-56.48 

17.10 

-26.33 

7 

19.52 

-57.07 

16.61 

-37.18 

9 

18.96 

-77.68 

15.70 

-33.38 

11 

15.45 

-81 .07 

14.95 

-50.38 

13 

15.58 

-77.76 

15.94 

-52.83 

15 

14.05 

-76.85 

14.27 

-56.08 

17 

12.19 

-77.49 

12.38 

-71 .38 

19 

11.50 

-68.77 

9.98 

-63 . 1 3 
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Table  H 2 - Describing  Function  Data  - Human  Subjects 
Frequency  Range:  0.05  Hz 


Subject 

- H 2 

Subject 

- H 3 

Freq.  (Hz) 
* 

122.6 

Mag.  (dB) 

Phase 
(Degrees ) 

Mag.  (dB) 

Phase 

(Degrees) 

1 

34.68 

-105.39 

45.18 

-59.07 

3 

30.31 

-64.19 

23.27 

-48.48 

5 

21.93 

-85.27 

19.37 

-57.19 

7 

23.72 

-108.01 

17.38 

-77.14 

9 

23.97 

-45.17 

19.15 

-63.91 

11 

23.82 

-85.42 

18.46 

-29.22 

13 

19.63 

-94.73 

16.56 

-89.22 

15 

18.46 

-71.19 

16.87 

-73.30 

17 

15.32 

-88.78 

14.46 

-69.00 

19 

12.84 

-108.24 

10.40 

-88.32 
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be  noted  that  the  data  presented  for  the  human  subjects  is  only  for  one 
run  each.  A mention  must  also  be  made  of  the  fact  that  even  though  the 
task  set  for  human  subjects  was  compensatory  in  nature,  it  was  similar, 
but  not  the  sane,  as  that  set  for  the  primates.  Whereas  the  primates 
were  punished  with  shock  for  their  off  target  performances,  the  human 
subjects  were  punished  for  these  errors  on  an  error-cost  basis  which 
subtracted  from  the  total  amount  of  money  ($3.00  maximum)  which  they 
could  earn.  Since  data  for  the  human  operators  were  only  recorded  for 
one  session,  no  statistical  analysis  was  performed.  The  data  for  these 
subjects  is  plotted  in  Figs.  2.14-2.17.  For  the  0.05  Hz  range,  both  for 
human  and  subhuman  operators,  ten  describing  function  estimates  were 
calculated  starting  at  the  fundamental  frequency  of  (1/245.2)  Hz  and 
repeated  for  9 frequency  points,  separated  in  frequency  domain  by 
(2/245.2)  Hz.  Similarly,  for  the  0.15  Hz  range,  the  describing  function 
estimates  were  computed  for  the  fundamental  frequency  of  (1/122.6)  Hz 
and  repeated  at  14  frequency  points  equally  spaced  at  (2/122.6)  Hz. 

Using  the  model  identification  scheme  discussed  earlier,  model 
parameters  for  all  subjects  were  calculated  and  are  reported  in  Table  P 
for  both  frequency  ranges,  .05  Hz  and  .15  Hz.  The  transfer  functions 
corresponding  to  these  model  parameters  are  plotted  in  Figs.  2.9-2.17, 
along  with  the  measured  describing  function. 

It  is  apparent  from  Tables  ‘12  and  M3  that  the  describing  functions 
measured  for  Butch  and  Hoppy  over  the  .15  Hz  range,  were  consistent  from 
day  to  day.  Also  Figs.  2.9  and  2.10  show  that  performances  of  Hoppy  and 
Butch  were  similar  over  the  .15  Hz  bandwidth.  Table  M4  shows  marked 
variations  in  Big  Boy's  performance  from  one  day  to  another.  This  in- 
consistenty  in  his  tracking  performance  is  also  confirmed  by  his  low 
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Fig.  2.15 
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correlation  index  p (id)  tabulated  in  Tables  R1  and  R2  respectively. 

“2 

The  p (id),  from  a control's  viewpoint,  is  a measure  of  the  linearity  in 

o 

the  operator's  response.  Alternatively,  [1  - p (id)]  is  a measure  of 

the  "remnant"  present  in  the  operator's  response.  A close  look  at  Tables 

2 

R1  and  R2  reveals  consistently  high  values  of  p for  Hoppy  and  Butch 

over  both  ranges  of  frequency  used.  The  values  typically  lie  between 

0.58  and  .96.  The  animals  showed  a marked  tendency  to  introduce  more 

remnant  at  higher  frequencies.  The  human  subjects'  responses  were  much 

more  linear  over  the  frequency  spectrum  used,  as  confirmed  by  high 
2 

p (id)  values. 

The  observed  inconsistency  in  Big  Boy's  response  can,  in  part,  be 

explained  in  terms  of  adaptation  and  motivation.  Again,  referring  to 

the  Remnant  column  in  Table  P it  may  be  noticed  that  the  reported 

values,  which  physically  represent  minimized  sum-squared-error  S,  are 

quite  consistent  for  all  subjects  except  Big  Boy,  The  least-squared 

2 

error  calculated  was  seemingly  low  and  inconsistent  with  p (w)  values. 
This  was  possible  due  to  a local  minima  observed  in  his  unpredictable 
response.  The  interanimal  consistency  as  seen  in  Hoppy 's  and  Butch's 
performance  could  be  attributed  to  their  ability  to  adapt  to  the 
tracking  task.  Adaptation  together  with  motivation  produced  the  same 
functional  relationships  regardless  of  the  species.  In  the  case  of  Big 
Boy,  motivation  seems  to  have  been  more  influential,  for,  even  though 
he  was  given  the  same  shock  treatment  every  day,  his  performance  was 
inconsistent  and  showed  large  variations.  In  the  .05  Hz  range  Hoppy 
showed  relatively  less  variation  than  Butch. 

The  data  for  human  subjects  plotted  in  Figs.  2.14-2.17  show  simil- 
arity of  trend  with  that  of  Rhesus  monkeys.  As  seen  from  Figs.  2.9-2.17 
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the  modeled  transfer  functions,  very  closely  match  the  measured  describ- 
ing function  both  for  human  and  subhuman  operators  except  Big  Boy. 

The  mathematized  and  quantitative  observations  discussed  above 
show  qualitative  agreement  with  preliminary  observations  made  on  the 
basis  of  digital,  Time-On-Target  scores  (TOT).  The  following  table 
shows  percent  time  on  target  for  all  subjects. 

COMPENSATORY  TASK 


Time-On-Target  (TOT) 


Frequencies 

Humans 

Monkeys 

0.05  Hz 

99.5 

94.7 

93.9 

99.0 

0.15  Hz 

99.3 

94.2 

98.5 

98.2 

90.8 

DISCUSSION 

1 . Input  Signal s 

Before  discussing  the  experimental  describing  functions  and  the  cal- 
culated transfer  functions  from  estimated  model  parameters,  it  is  appro- 
priate to  examine  the  nature  of  the  time  signals  circulating  in  the 
operator/machine  loop.  The  overall  system  inputs  were  derived  from  a 
Gaussian  noise  generator  and  the  operator  attempted  to  force  the  system 
output  to  equal  the  system  input  by  operating  on  the  stick  controller. 
Depending  on  the  set  performance  criteria,  operator's  output  may  or  may 
not  have  been  Gaussian.  To  have  a better  insight  into  this  aspect  of 
the  problem,  the  width  of  the  window  will  have  to  be  taken  into  account 
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(see  Window  Width).  However,  since  the  input  signal  x(t)  was  an  extern- 
ally generated  random  noise,  it  would  be  uncorrelated  with  any  noise 
generated  within  the  system  and  present  at  the  output.  This  precisely 

9 

means  that  the  correlation  index  p ' (w)  given  by 

o k...  U)'2 

2.18 


p2  (uj)  = 


! xy 

<J>  (oj)  <l>  (w) 

xx  ' ' yy 


would  be  a true  measure  of  the  nonlinearity  in  the  operator  response. 

Apart  from  being  Gaussian-random,  the  noise  bandwidth  was  set  to  a 
maximum  cutoff  frequency  of  .15  Hz.  This  extremely  low  frequency  range 
was  primarily  used  for  two  reasons:  first,  to  train  the  animals  on  a 
simple  task  and  then  gradually  make  the  task  more  difficult.  The 
second  reason  was  to  provide  a data  base  in  a frequency  range  that  has 
hitherto  been  uninvestigated.  Further,  it  was  observed  that  at  higher 
frequency  ranges,  the  animals  would  not  track  despite  the  shock  treatment. 


2.  Window-width 

The  error  signal  e(t)  was  generated  by  taking  the  difference  of 
y(t)  and  x(t).  This  would  have  been  mathematically  appropriate  if  the 
subjects  were  required  to  align  the  cursor  with  a central  line  on  the  CRT 
display.  In  such  a situation  a leftward  movement  of  the  cursor  would 
have  required  rightward  compensation  and  vice-versa.  Since  this  task 
would  have  been  extremely  difficult,  the  subjects  were  required  to  align 
the  cursor  with  a line  of  finite  width  (one  inch).  From  a control's 
standpoint  this  changed  the  nature  of  the  task.  Quantitatively,  the 
error  signal  e(t)  could  always  have  been  non-zero  even  though  the  sub- 
ject was  tracking  99%  of  the  time.  In  principle,  the  animal  was  not 
given  a true  indication  of  the  error  signal.  Consideration  must  also 
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be  given  to  the  fact  that  the  animal  could  have  just  been  sitting  and 
wiggling  the  stick  generating  error  even  though  his  TOT  scores  indicated 
zero  or  no  error.  Conclusively  then,  simplification  in  task  resulted 
in  ambiguities  in  operator  performance.  This  is  confirmed  by  Big  Boy's 
performance  results.  His  TOT  scores  were  appreciably  high  (90.0%)  even 
though  the  correlation  index  reflected  a very  poor  performance  (.30). 

3.  Human  vs.  Subhuman  Operators. 

As  part  of  the  research,  data  were  collected  on  human  and  subhuman 
subjects  both  in  compensatory  and  pursuit  tracking  situations.  However, 
in  the  case  of  pursuit,  since  only  discrete  input  frequencies  (.122  Hz 
and  .163  Hz)  were  used,  describing  function  data  was  not  calculated. 
Therefore,  comparisons  between  the  two  species  will  be  made  solely  on  the 
basis  of  results  obtained  from  the  compensatory  tracking  task.  In  this 
context,  it  is  of  interest  that  studies  reported  by  Bachman  et  al  . 

(1976)  focused  on  the  issue  of  man-monkey  comparisons  in  a tracking 
situation.  They  used  a primate-training  task  quite  similar  to  that  dis- 
covered in  our  laboratory  a few  years  ago.  Bachman  et  al . have  re- 
ported results,  similar  to  those  observed  in  this  investigation,  which 
lend  credibility  to  monkey-man  extrapolation.  Bachman  et  al.  observed 
very  low  "remnant"  in  monkey  response  which  is  essentially  what  is  re- 
ported in  Tables  R1  and  R2. 

Quantitative  evaluation  of  man-monkey  performance  revealed  that  for 
both  the  human  and  subhuman  operators  the  transportation  lag  was  within 
.08-. 19  seconds.  This  is  comparable  with  transportation  delay  reported 
in  literature  for  human  operators  in  similar  tracking  situations,  Shinners 
(1974),  McRuer  et  al . (1974)  and  Bachman  et  al . (1976). 


58 


2 

Table  R 1 - Correlation  Index  - p (w) 
Frequency  Bandwidth  - 0.15  Hz 


donkeys 


Human  Subjects 


Freq.  (Hz) 


★ 

122.58 

Big  Boy 

Butch 

Hoppy 

H 2 

H 3 

1 

.23 

.92 

.89 

.99 

.99 

3 

.28 

.94 

.94 

.99 

oo 

as 

5 

.40 

.94 

.93 

.98 

co 

Os 

7 

.27 

.95 

.93 

.99 

.98 

9 

.45 

.96 

.94 

.99 

.98 

11 

.32 

.95 

.93 

.99 

.98 

13 

.35 

.95 

.89 

.98 

.98 

15 

.24 

.95 

.92 

.99 

.90 

17 

.23 

.95 

.90 

.98 

.98 

19 

.19 

.94 

.93 

.99 

.99 

21 

.25 

.93 

.90 

.98 

.97 

23 

.22 

.89 

.83 

.98 

.96 

25 

.20 

.86 

.80 

.92 

.85 

27 

.16 

.80 

.76 

.88 

.82 

29 

.14 

.79 

CO 

.89 

.88 

r 
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Table  R 2 - Correlation  Index  - p (a>) 
Frequency  Bandwidth  - 0.05  Hz 


Freq.  (Hz) 
★ 

245.2 

Monkeys 

Human 

Subjects 

Butch 

Hoppy 

H 2 

H 3 

i 

.77 

.91 

.98 

.99 

3 

.78 

.94 

.98 

.99 

5 

.77 

.94 

.99 

.99 

7 

.77 

.82 

.98 

.98 

9 

.78 

.82 

.98 

.99 

11 

.76 

.82 

.98 

.99 

13 

.79 

.87 

.97 

.97 

15 

.76 

.81 

.97 

.97 

17 

.81 

.66 

.91 

.97 

19 

.87 

.58 

.91 

.95 

A 
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The  results  obtained  in  this  investigative  study  along  with 
those  of  Bachman  et  al . offer  strong  convincing  evidence  for  the  con- 
tention that  both  these  species  show  similar  performances.  Therefore, 
the  fundamental  differences  between  man  and  monkey  reported  Ly  some 
do  not  match  empirical  results  obtained  in  this  study.  Further  exam- 
ination of  this  will  require  further  experimentation  and  more  complex 
tracking  tasks. 

4.  Comparison  with  Data  in  the  Literature. 

To  this  author's  knowledge,  this  kind  of  research  usinq  primate 
operators  has  only  been  reported  by  Bachman  et  al . (1976).  This  re- 
search was  partly  conducted  to  add  to  the  rather  meager  data  base  of 
primate  tracking  experiments.  The  table  below  shows  "best  fit"  para- 
meters reported  by  Bachman  et  al . 


Table  L - Data  from  Bachman  et  al . 


Parameter 

Subj.  H 1 

Subj.  H 2 

Subj . M 1 

Subj.  M 2 

K,  Gain 

20.0 

11.0 

7.5 

8.0 

T^,  Lead 

4.2 

1 .4 

4.2 

1 .6 

V Lag 

41.9 

31  .4 

20.9 

20.9 

Tj,  Lag 

1 .6 

0.3 

2.1 

0.5 

t.  Time 
Delay 

0.14 

0.21 

0.14 

0.19 

Even  though  the  results  of  the  reported  study  (see  Table  P)  reflect  re- 
markable similarity  with  the  observations  made  by  the  above  authors, 
considerable  care  must  be  exercised  in  interpreting  the  results.  Re- 


sults of  Table  P show  less  variation  than  those  in  Table  L above.  This 
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Table  P - Model  Parameters 


Freq . 
Range 
(Hz) 

Subject 

Name 

T 

(sec) 

T1 

(sec ) 

T2 

(sec) 

T3 

(sec) 

K 

0.05 

H 2 

.16 

18.01 

28.99 

3.99 

13.97 

H 3 

.16 

18.00 

28.99 

3.99 

18.77 

Butch 

.09 

18.01 

28.99 

4.59 

13.00 

Hoppy 

.09 

18.01 

28.88 

3.99 

13.00 

0.15 

H 2 

.09 

18.00 

28.99 

3.91 

13.70 

H 3 

.17 

18.00 

27.10 

3.98 

18.98 

Big  Boy 

.19 

18.00 

28.98 

3.99 

13.00 

Butch 

.10 

16.50 

24.36 

4.5 

11.99 

Hoppy 

.08 

18.00 

25.00 

3.99 

18.99 

Remnant 


418. 

288.3 

398.4 
365. 0 

22.9 

276.7 

187.9 

25.1 

336.6 


r-o 
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is  due  to  the  least  squared  error  optimization  used  in  this  study.  The 
parameteric  variation,  to  some  extent,  is  reflected  in  the  Remnant  col- 
umn of  Table  L.  The  above  authors  used  the  BMD  02T  computer  program  for 
estimating  describing  function;  Fig.  2.7  shows  the  superiority  of  com- 
puter algorithm  used  in  this  study  over  BMD  02T.  Also  their  training 
procedure  was  similar  to  the  one  used  in  this  study,  but  the  paper  fails 
to  bring  out  the  ramifications  of  the  window  size.  Finally,  the  optimi- 
zation algorithm  of  Section  E has  a definite  superiority  over  visual 
"interactive-fit"  used  in  the  above  reference. 
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